On the large magnetic anisotropy of Fe 2 P 
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We present an investigation on the large magnetic anisotropy of Fe2P, based on A b Initio density- 
functional theory calculations, with a full potential linear muffin-tin orbital (FP-LMTO) basis. We 
obtain an uniaxial magnetic anisotropy energy (MAE) of 664 fieV/i.u., which is in decent agreement 
with experimental observations. Based on a band structure analysis the microscopical origin of the 
large magnetic anisotropy is explained. We also show that by straining the crystal structure, the 
MAE can be enhanced further. 



INTRODUCTION 

Magnetic materials with large magnetic anisotropy 
(MAE) have been used in many applications. Until re- 
cently the dominating class of materials used for such 
applications, apart from hard ferrites, were rare-earth 
based magnets, e.g. Nd2Fei4B in which the large 
magnetic anisotropy energy is provided by the rare-earth 
atoms, and the large saturation moment (M s ) at finite 
temperature is due to the Fe atoms. It has however been 
pointed out that other permanent magnets should be in- 
vestigated, from an application point of view, since a 
general access to rare-earth elements is far from guaran- 
teed 0. Among such materials the iron-phosphide Fe2P 
stands out as a particularly interesting material, due its 
known large MAE and sufficient large value of satura- 
tion moment . In addition Fe2P is composed of cheap 
and widely available elements. 

Transition metal pnictidcs and chalcogenides are not 
the only large MAE materials that have been investi- 
gated so far. For instance, a large saturation moment 
was suggested in a nano-laminate of a 3d metal (Fe) and 
a rare earth metal (Gd) in ReflH. Also, a large mag- 
netic anisotropy in FePt was verified experimentally as 
well as from first principles theory Furthermore, 
the predicted large MAE of a tetragonally strained FeCo- 
alloy Q was verified experimentally (Io| . 

In this report we focus on the MAE of Fe2P, since it is 
large and hitherto unexplained. In order to find a micro- 
scopic description of the large MAE of this compound we 
have performed first principles calculations, using a rel- 
ativistic formulation of the Kohn-Sham equation. There 
arc several experimental studies concerning the ordering 
temperature, MAE, saturation moment, hyperfinc field 
and isomer-shift of Fe2P (llT - flil ]. In addition, a theoret- 
ical analysis was made earlier by Wohlfarth [ljj [lj| , and 
a subsequent theoretical work addressed the magnetism 
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of Fe2P using electronic structure calculations 
An excellent overview of the magnetic properties of FC2P 
and similar transition metal pnictides and chalcogenides 
can be found in Refflii 

In its crystal structure Fe2P (Hexagonal C22 with space 



group P62m, No. 189) [20|,[2l| has two Fe sites. The Fe-I 



site has a magnetic moment close to 1 /is/atom, and the 
Fe-II site has a magnetic moment close to 2 /xs/atom (In 
Fig. [T]the structure of Fe2P is shown). The material has 
an observed MAE of 500 /ueV/f.u. [HI , with the crystal- 
lographic c-axis being the magnetization easy axis. The 
microscopic mechanism for this large MAE is however 
unknown. It is the motivation of the present study to 
find this mechanism, and to suggest Fe2P based alloys 
with enhanced values of MAE. 




Fig. 1: (color online) Fe2P crystal structure. Fe-I tetrahedral- 
sites (red), Fe-II pyramidal-sites (green) and P atoms (yel- 
low). 



COMPUTATIONAL DETAILS 

We used a full potential linear muffin-tin (FP-LMTO) 
method (22[, with a relativistic formulation, where the 
spin-orbit coupling is included inside the muffin-tins 
spheres, at the variational step. The muffin-tins are 
centered on atomic sites, and a combination of radial 
functions multiplied by spherical harmonics are used to 
expand the electronic density and potential inside each 
muffin-tin. In the interstitial region a combination of 
Hankel and Neumann functions are used as basis func- 
tions. We calculated the MAE by the magnetic force the- 
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orem 23| . A scalar relativistic calculation, with common 
symmetry of both magnetization directions, is first per- 
formed to obtain a self-consistent potential. With this 
potential fixed one performs a fully-relativistic calcula- 
tion for each magnetization orientations (hi). The Emab 
(magnetic anisotropy energy) is then calculated from the 
expression 



Unfortunately most dopings on the P site, e 
or B, results in a reduced c/a ratio [l2, 2(| 2f 



with Si 



Emae = ^ e *("2, k) 



i(hi, k). 



(1) 



Here ej(n, k) is the Kohn-Sham eigenvalues evaluated for 
each magnetization orientation. In Eq.l, i labels the oc- 
cupied states and k the k-points in the Brillouin zone, 
whereas hi = 0001 (c-axis) and n2 = 1000. 

The modified tetrahedron method (MTM) [H] was 
used for integration in the Brillouin zone (BZ), 36002 k- 
points were used (for a calculation with 50626 k-points 
the MAE deviates only 0.1 %), in irreducible part of BZ, 
to guarantee the convergence. 



RESULTS 

Our main result is shown in Fig. [2J where we display 
the calculated MAE as a function of the strain. Note that 
we have used the room temperature lattice constant as 
reference level, having room temperature applications in 
mind. It is clear from the figure that theory reproduces 
the observed easy axis (0001) and that the calculated 
MAE is of the same order of magnitude as the experi- 
mental one. The theory overestimates the value of the 
MAE by 32 %, when a comparison is made between the 
experimental low temperature value (red square, taken 
at 4 K) and the lattice parameters corresponding to this 
temperature (which in Fig. 2 corresponds to a - 0.55 %). 
The strain state of the low temperature lattice constant 
is calculated from the work of Fujii et. 
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ported measurements for the thermal expansion of the 
FC2P in the temperature range of 60 to 550 K. Using 
this thermal expansion data the 4 K lattice parameters 
were estimated. The agreement between calculated and 
measured MAE obtained here is typical, when compared 
to other calculations [fj H, 25 1, and primarily reflects the 
extremely delicate nature of the MAE in general. A very 
important result shown in Fig. [21 however, is that with an 
applied strain to the lattice it is possible to influence the 
MAE quite substantially. In these calculations, we ap- 
plied strain, keeping a constant volume. Increasing the 
c value by 1% affects the MAE to be enhanced by ~ 15 
%, which is the maximum value of the calculated MAE. 
For strain values higher than 1%, Fig. [2] shows that the 
MAE decreases almost linearly. Reducing the c-axis is 
not favorable for the MAE, here we note a monotoni- 
cally decreasing trend, with a minimum observed at -8%. 



The magnetic moment of Fe2P changes much less with 
strain when compared to the MAE, as the inset of Fig. [21 
shows. At 1% strain the calculated moment is 3.03 
//s/f.u., which is -0.4 % of the zero strained system. The 
magnetic moment increases however slightly for negative 
values of strain. Overall our calculated moments agree 
well with the observed number of 2.94 /xs/f.u 
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The Fig. [2] shows data for a volume conserving strain 
(except the red square which corresponds to a volume 
which is minutely smaller than the volume used for the 
other data points). For comparison we also show in Fig. [3] 
the MAE as function of strain in two non- volume conserv- 
ing regimes. 

This involves strain of the c-axis while a and b are fixed 
as well as strain of the ab-axis while c axis is fixed. For 
c-axis (ab-axis) strain a maximum MAE value of 800 /icV 
(807 fieV) is obtained for 2% (-1%) strain. Experimen- 
tal doping on the P site with B results in a reduced c/a 
axis (negative strain) and a reduced volume [l2. 2fl 26j] . 
in this case a predicted MAE with doping is best eval- 
uated by inspection of the non- volume conserving curve 
in Fig. [3l Doping with Si keeps the volume essentially 
constant while reducing the c/a ratio, and for this dop- 
ing element the predicted MAE is best evaluated from 
the volume conserving curve, which shows a decreased 
MAE. The MAE was actually calculated using the struc- 
tural cell parameters for 10 % Si doping [l^|, and 593 /xeV 
MAE was obtained, showing the expected reduction. For 
10% Si doping the experimental Curie temperature (T c ) 
is 370 K (an increase of 70 % when compared with the un- 
doped case) , which is promising for stabilizing a material 
which at room temperature has a large MAE. Further in- 
vestigations are necessary to consider the chemical effect 
impact on the MAE. 

Both these predicted changes of the MAE with doping 
rely on the applicability of the rigid band approximation. 
They have not been evaluated experimentally, and a ver- 
ification or refutal of this prediction would be interesting. 

The strong variation of the MAE with strain in Fig. [2] 
is interesting and requires further analysis. On a simple 
model level, the uniaxial MAE can be treated in second- 
order perturbation theory, and computed as the differ- 
ence of the second order correction to the energy (E^ s [h]) 
between two magnetization directions as in Eq. [2J with 
a sum over the sites q, over occupied (s) and unoccupied 
(s') spin characters, respectively; [27 1 



AE so = AE f = EW("a) - Efih!)}, (2) 



qss' 



where the energy correction is given by, 
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Fig. 2: (color online) Calculated MAE as function of the 
strain. Positive values of the strain correspond to an increased 
c/a ratio (at constant volume). Positive values of the MAE 
correspond to the c-axis being the easy axis. The inset shows 
the calculated magnetic moment as a function of strain. The 
square (red) is the MAE value for the estimated 4 K lattice 
parameters. 



basis functions \qlms) are characterized by the site q, 
azimuthal (1), magnetic (m) and spin (s) quantum num- 
bers and £ki (fkj) are the electronic eigenvalues for the 
occupied (unoccupied) states. The hybridization is con- 
sidered in the band character nkis.gm,g'm' , which allows 
mixing of basis functions on different sites. 
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Fig. 3: (color online) Calculated MAE as function of the 
strain of the c-axis with a and b fixed (red squares) and strain 
on the ab-axis with c fixed (black circles). 



(Ti) — ^ ^ ^ ^kzs,qm,g / m / ^kjs / ,q'm" ^qra"' X 
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£kj — £ki 



(3) 



In Eq.[3]thcrc is a sum over k points in the Brillouin zone, 
i and j label the occupied and unoccupied states, s and 
s' run over the spin character of the states and m, m', 
m", m'" run over the magnetic quantum numbers. The 



Fig. 4: (color online) Band structure of Fe2P calculated with- 
out spin-orbit coupling (the band thickness represents the 
weight of Fe-II, 1=2 and m=-2 state). Red (solid) is the spin 
up character and blue (striped) is the spin down. The Fermi 
level (Ef) is placed at E=0. 



Since the electronic eigenvalues (eki) appear in Eq. [3l 
it is relevant to inspect the band structure along the high 
symmetry lines of the hexagonal lattice. Hence, we show 
the calculated energy bands in Fig. @J where the thick- 
ness of the bands represents the weight of the Fe-II, 1=2 
and m=-2 state. The bands that contribute the most to 
the MAE are highlighted by the arrows 1 to 5. At |(K- 
M), arrow 1, the occupied and unoccupied bands have 
mainly character from 1=2 and m=±2 quantum num- 
bers, and these bands interact through the l z s z term of 
the spin-orbit coupling (SOC) Hamiltonian, resulting in 
a large negative contribution. A similar mechanism is 
observed at |(r-A), arrow 4. This is illustrated further 
using a calculation with spin-orbit coupling included, in 
Fig. [5] a), in a region zoomed in around the Fermi level, 
see the highlighted areas. One can see the splitting of the 
bands when the spin quantization axis is along the 0001 
crystallographic direction. The splitting of these bands 
is not observed for the 100 axis, as observed in Fig. [5] 
b). At arrows 2, 3 and 5 the occupied and unoccupied 
bands have different m quantum numbers, m=±2 and 
m=0 (m=±l) for the occupied (unoccupied) bands, the 
SOC interaction via the terms : (+s_+Ls+, give large 
positive contributions. As a general rule if the occupied 
and unoccupied bands have the same (different) m quan- 
tum number, the magnetization is favored to lie along 
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Fig. 5: (color online) Zero strain band structure of Fe2P 
calculated with spin-orbit coupling, magnetization along the 
a) 001 and b) 100 crystallographic axis. The Fermi level (Ef) 
is placed at E=0. 



the 001 (100) axis j28|. 

To understand the enhancement of the MAE under 
uniaxial strain, the bandstructure was calculated for the 
1 % volume conserving strain (Fig [5] a). The overall 
bandstructure is not changed, due to the strain, only 
a rigid shift is observed. Without SOC the flat bands 
around the A point lies on the Fermi level. As pointed 
out before these bands are split due to SOC for a mag- 
netization along the 001 axis. One of the bands become 
fully occupied and the other fully unoccupied, given that 
these bands are very close in energy a large contribu- 
tion to the MAE is expected. Indeed this is confirmed 
in Fig |6] b) (upper panel) where the MAE for each k- 
point (EjvfA_E[k]) along the high symmetry lines of the 
hexagonal cell is showed, for the zero and 1% (volume 
conserving) strained system. For the zero strain one can 
see a sharp negative peak for the MAE (negative val- 
ues favors the 001 axis) around arrows 1 and 4. For 1% 
strain the peak at arrow 4 gets broadened explaining the 
increasing of the MAE. 
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Fig. 6: (color online) a) Band structure of Fe2P calculated 
with 1 % of strain on the c-axis and b) MAE along the high 
symmetry directions. The Fermi Level (Ef) is placed at E=0. 



CONCLUSION 

In conclusion we have studied the uniaxial MAE of 
Fe2P by Ah Initio calculations. Our theory repro- 
duces the observed 0001-easy axis with a MAE of 664 
/ieV/f.u. This should be compared to an experimen- 
tal value of - 500 ^eV/f.u.(2.32 MJ/m 3 ) This is an 
acceptable agreement between theory and experiment, 
when having in mind the extremely delicate nature of 
the MAE and the typical low energy differences associ- 
ated with it. The size of the MAE of Fe 2 P should be 
compared to other hard magnetic materials like FePt 1.2 
meV/f.u. (6.6 MJ/m 3 ) Q and Nd 2 Fe 12 B 6.7 meV/f.u. 
(4.9 MJ/m 3 ) [2^]. As to the magnetic moments our cal- 
culations give value of 3.04 /Ltg/f.u., which agrees well 
with the experimental moment of 2.94 /is/f.u. 

We have analyzed the origin of the MAE by a detailed 
band and k-point resolved property, and show that the 
positions of the different energy bands around the Fermi 
level critically determine that MAE. Since these bands 
can be moved up or down in energy with an applied 
strain, it seems that this is an important avenue with 
which to influence the MAE of Fe 2 P and Fe 2 P-based al- 
loys. Consequently we find from our MAE calculations 
that it is possible to influence the MAE quite strongly 
with an applied strain, both in a volume conserving and 
a non-volume conserving mode. This opens up possibili- 
ties to use Fe 2 P and alloys of this material, as a platform 
for searching for new permanent magnetic materials that 
don't contain rare-earth elements. 
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